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1 Introduction and framework 

oo: 

G^ ■ 1.1 Overview 
O; 

CO I The well-known Fortmanteaw r/ieorem of measure theory (see e.g. Op. 15 ff.]) is a powerful statement, 

■ providing several necessary and sufficient conditions in order for a sequence of probability measures on 
I a metric space to converge weakly towards some limit. The term 'portmanteau' indicates that these 

04 ■ conditions have a priori different natures, in such a way that they appear as artificially packed together 

' at first readingH 
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The aim of this paper is to prove and apply a new portmanteau inequality, involving vectors of random 
variables that are functionals of a Poisson measure defined on a general space. This estimate - which 



H . is formally stated in formula (|2.9I) below - is expressed in terms of Malliavin operators, and basically 

allows one to measure the distance between the laws of a general random element and of a random vector 
whose components are in part Gaussian and in part Poisson random variables. As we shall abundantly 
illustrate in the sequel, the inequality ()2.9p is a genuine 'portmanteau statement' - in the sense that it can 
be used to directly deduce a number of new results about the convergence of random variables defined 
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on a Poisson space, as well as to recover known ones. These results span a wide spectrum of asymptotic 
behaviors that are dealt with for the first time in a completely unified way. Apart from Malliavin calculus 
(that we apply in a form analogous to the one developed by Nualart and Vives in [34]), our techniques 
involve the use of the Chen-Stein method (see e.g. [3]), and provide a substantial refinement of several 
recent contributions concerning Central Limit Theorems (CLTs) and Poisson approximation results on 
the Poisson space (see 27, 28, 35l [37l HH [5TJ [57]). One of our main technical tools is an interpolation 
technique used in [3] for proving multidimensional Poisson results. See e.g. [32l [33] for a discussion of 
the use of Stein-Malliavin techniques on a Gaussian space. 

As the title indicates, the two new main theoretical applications developed in the sequel are the following: 

— Mixed limits: Our results allow to deduce quantitative limit theorems (that is, limit theorems 
with explicit information on the rate of convergence), where the target distribution is a multi- 
dimensional combination of independent Gaussian and Poisson components. This new class of 
approximation results is described in Section 12.11 They will be applied both to characterize the 
asymptotic independence of general ?7-statistics (see Section 12. 3p , and to subgraph counting in 
stochastic geometry (see Section 12.41) . By virtue of an approximation argument borrowed from 
[16], part of the results discussed in Section [2731 extends to de-poissonized [/-statistics. 

- Multi- dimensional Poisson convergence: A particular choice of parameters in our main estimates 
allows one to deduce multidimensional Poisson approximation results, having moreover a stable 
nature - in the classic sense of [2] [S^. This generalizes the one-dimensional findings of [35]. 
See Section 12.2.11 and Section I2.2.2[ respectively, for general statements and for applications to 
sequences of multiple Wiener-Ito integrals, as well as for several comparisons with the CLTs 
established in [ST] [44] . Characterizing the convergence in distribution of random variables having 
a chaotic nature (both in a classic and a free setting) has recently become a relevant direction 
of research (see e.g. [33j for an overview of the many available results in a Gaussian setting, or 
[HI [5S] for several free counterpartf[l) , and our analysis provides substantial new contributions in 
the case of random variables belonging to the Poisson Wiener chaos. One should also note that 
Poisson approximation results based on Malliavin operators have found a number of applications 
in stochastic geometry, see [59] . 

The basic intuition underlying our approach is the following: in order to properly understand the 
connections between Poisson approximations and CLTs in the context of random point measures, it 
is very much instructive to study probabilistic models where Poisson and Gaussian random structures 
emerge simultaneously in the limit. The present paper demonstrates how Portmanteau inequalities 
provide the correct tool for accomplishing this task in a fully multidimensional setting. 

We will illustrate our findings by completely developing an application to random geometric graphs, as 
described in Section [2741 and Section [4] In particular, two results will be achieved: 

(i) a new bound for the multidimensional Poisson approximation of subgraph-counting statistics; 

(ii) a proof of a new mixed limit theorem involving the joint convergence of vectors of subgraph- 
counting statistics exhibiting both a Poisson and a Gaussian behavior. 

Our results extend several findings in the field ~ see [5l [23l [45] . 

complete list of the papers related to this topic can be retrieved from the constantly updated webpage 
http: //www. iecn.u-nancy . f r/~nourdin/steinmalliavin.htm 
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Remark 1.1 Due to the use of the Chen-Stein method, one of the main technical difficulties in develop- 
ing our results has been the choice of a "good" version of a discrete multivariate Taylor-type expansion 
for functions defined on the set Z^j. = {0, 1, 2, ...}'', d > 2. The formula that best fits our approach ap- 
pears in Lemma |3. II it provides a representation of the remainder as a double sum, where diagonal and 
non-diagonal terms play asymmetric roles. Our analysis implicitly shows that such a formula virtually 
encodes all the combinatorial subtleties involved in the derivation of Poisson approximation results on 
the Poisson space. See Barbour for several applications of univariate discrete Taylor formulae to the 
computation of factorial moments and cumulants. 

The remainder of the paper is organized as follows. The next subsection contains a formal description 
of our framework: it is mostly standard material, so that someone already familiar with the notation of 
[2711281 [571 144) can skip it at first reading. Section [5] contains a detailed discussion of the main theoretical 
results of the paper, as well as of the applications. Section [3| is devoted to the proofs of our general 
theorems, whereas Section |4| contains the proofs of our results about random graphs. An Appendix 
contains basic notions about Malliavin operators and contractions. 

1.2 Framework 

In what follows, we shall denote by {Z, 3f, a measure space such that Z is a Borel space, 2f is the 
associated Borel tr-field, and is a cr-finite Borel measure with no atoms. We write = {B G Z : 
IJl{B) < oo}. The notation rj ~ {viB) ■ B e i^} is used to indicate a Poisson measure on {Z,Z) 
with control (or intensity) ^. This means that 77 is a collection of random variables defined on some 
probability space (^2, P), indexed by the elements of 2ff_i and such that: (i) for every B,C £ such 
that J5 nC — 0, the random variables ri{B) and 77(C) are independent; (ii) for every B g v{B) has 
a Poisson distribution with mean fJ.{B). We shall also write 

^%B)=^^{B)-^I{B), B e 

and 77 = {viB) ■ B e A random measure verifying property (i) is usually called "completely 

random" or "independently scattered" (see e.g. [42j [55] for a general introduction to these concepts, 
and for a discussion of any unexplained definition or result). 

Remark 1.2 (The probability space) (i) In view of the assumptions on the space {Z, 3f, n), and 

to simplify the discussion, we will assume throughout the paper that (r2,^,P) and 77 are such 
that 

^ " 1^ = X^f^^,,?! e Nu {oo},zj e zj , 

where Sz denotes the Dirac mass at z, and 77 is defined as the canonical mapping 

{uj,B)^ri{B){uj)^Lo{B), B e JT^, a; G 17. 

Also, the cr-field ^ will be always supposed to be the cr-field generated by 77, and we will write 
I/^(P) = L^(ri,,^,P). Note that the fact that ^ is non-atomic implies that, for every x G Z, 
¥{ri{x} or 1} = 1 . 

(ii) As usual, by a slight abuse of notation, we shall often write a; <E 77 in order to indicate that the 
point X £ Z IS charged by the random measure 77(-). 
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Throughout the paper, for p E [1, oo), the symbol L^ifJ.) is shorthand for L'''{Z, 3f, /i). For an integer q > 
2, we shah write := ^"^9, /i*), whereas if (m') stands for the subspace of LP(/Lt'J) composed 

of functions that are /i'^-almost everywhere symmetric. Also, we adopt the convention = LP{^) = 

LP{fi^) — LP{ii^) and use the following standard notation: for every q>l and every /, g G L'^{fi'^), 

{f,9)L^f,l) ^ f{zi,...,Zq)g{zi,...,Zq)fJ,'^{dzi,...,dZq), |1/|1l2(^,,) = 

For every / e L?{ijfl), we denote by / the canonical symmetrization of /, that is, 

where a runs over the q\ permutations of the set {!,..., q}. Note that ||/||L2(^g-) < ||/||L2(^g-) (to see 
this, use for instance the triangular inequality) . 

Definition 1.1 For every deterministic function h e L^in), we write 

h{h)=f^{h) = I h{z)fi{dz) 
Jz 

to indicate the Wiener-Ito integral of h with respect to fj. For every q > 2 and every / e 

we denote by Iq{f) the multiple Wiener-Ito integral, of order q, of / with respect to 77. We also set 

Iq{f ) — Iqif), for every / £ L?'{p'^) (not necessarily symmetric), and /o(fo) = b for every real constant h. 

The reader is referred for instance to ^42^ Chapter 5] or [49] for a complete discussion of multiple 
Wiener-Ito integrals and their properties (including the forthcoming Proposition 11.11 and Proposition 

Proposition 1.1 The following equalities hold for every q,m > 1, every f G L^(/i'^) and every g G 
1. E[/,(/)] =0, 

2- Wq{f)Irn{9)] = q^f,9)L^t,i)'^{q=ni} (isometric property). 

The Hilbert space composed of the random variables of the form Iq{f), where q> I and / G L'^d.i'^), is 
called the qth Wiener chaos associated with the Poisson measure 77. The following well-known chaotic 
representation property is an essential feature of Poisson random measures. Recall that ^ is assumed 
to be generated by 77. 

Proposition 1.2 (Wiener-Ito chaotic decomposition) Every random variable F G L^(P) admits 
a ( unique ) chaotic decomposition of the type 

00 

F = nF]+Y,h{fi). (1.1) 

1=1 

where the series converges in L'^(P) and, for each i > 1, the kernel fi is an element of L^{p''). 



4 



Remark 1.3 (About Malliavin calculus) For the rest of the paper, we shall use definitions and 
results related to Malliavin-type operators defined on the space of functionals of the Poisson measure rj. 
Our formalism is the same as in Nualart and Vives in [34,. In particular, we shall denote by 

D, S, L and L~^, 

respectively, the Malliavin derivative, the divergence operator, the Ornstein- Uhlenbeck generator and 
its pseudo-inverse. The domains of D, S and L are denoted by domZ^, dom(5 and domL. The domain 
of L^^ is given by the subclass of i^(P) composed of centered random variables. For the convenience 
of the reader we have collected some crucial definitions and results in Section lA.ll of the Appendix. 
Here, we just recall that, since the underlying probability space fi is assumed to be the collection 
of discrete measures described in Remark 11.21 then one can meaningfully define the random variable 
uj I— > Fz{uj) = F{lu + 6z), w S ri, for every given random variable F and every z d Z, where is the 
Dirac mass at z. One can therefore prove the following neat representation of D as a difference operator 
is in order: for each F G dom_D, 

D^F = Fz - F, a.e.-n{dz). (1.2) 

Observe that the notation Fz{uj) = F(uj + 6z) extends canonically to multivariate random elements. A 
complete proof of this point can be found in ^4] . 

The next statement contains an important product formula for Poisson multiple integrals (see e.g. |42j for 
a proof). Note that the statement involves contraction operators of the type ★J,: the reader is referred to 
Appendix I A . 2 1 for the definition of these operators, as well as for a discussion of some relevant properties. 

Proposition 1.3 (Product formula) Let f €E Ll{iJ') and g G Ll{^'^), p,q>l, and suppose moreover 
that f ^J, g G L^(/,t^+^^''^') for every r — 1, . . . ,p A q and I — 1, ... ,r such that I ^ r. Then, 

ip{f)i,{9) = E^! ( ^ ) ( ' ) E ) (Mg) . (1-3) 

r=0 ^ ^ ^ ^ i=0 ^ ^ 

with the tilde ^ indicating a symmetrization, that is, 

f *l- 9{xi, . . . , Xp+q-r-l) = ^ q J. ^)!'^ ^'■^'^(1)' ■ ' ■ ' •'^<^iP+Q~r~l)), 

where a runs over all {p + q — r — ly. permutations of the set {1, . . . ,p + q — r — I}. 

Assumption 1.1 (Technical assumptions on kernels) In the sequel, whenever we consider a ran- 
dom vector of the type 

(/,,(/i),...,/,,(/d)), where d>l, > 1, /. G l2(/,?0, 

we will implicitly assume that the following properties (l)-(3) are verified. 

(1) For every i = I, d and every r = 1, qi, the kernel fi *qi~^ fi is an element of L^(/i''). 

(2) For every i such that f/i > 2, every contraction of the type (zi, Z2gi-r-i) -Z2<}i-r-i) 
is well-defined and finite for every r = l,...,qi, every I = l,...,r and every (zi, Z2q--r-i) G 
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(3) For every i,j — 1, ...,d such that max(gi, qj) > 1, for every k ~ |gi — | V 1, + — 2 and 
every (r, I) verifying k~qi + qj~2~r~l, 



li{dz) < CO, 



where, for every fixed z G Z, the symbol fi{z, •) denotes the mapping (zi, Zq^i) i— fi{z, zi, Zq-i). 

Remark 1.4 According to 031 Lemma 2.9 and Remark 2.10], Point (1) in Assumption 11.11 imphes that 
the following properties (a)-(c) are verified: 

(a) for every 1 < i < j < k, for every r — 1, qi A qj and every I = 1, r, the contraction fi ik-J. fj is 
a well-defined element of L^(^'''+'?J^'"^'); 

(b) for every I < i < j < k and every r = 1, qi, fi -k'^ fj is an element of 

(c) for every i = 1, k, for every r — 1, q^, and every I = 1, r A {qi — 1), the kernel ^J, fi is a 
well-defined element of 

In particular, every random vector verifying Assumption 1 1.1 1 is such that Iq-{fi)'^ G L^iP) for every 
i = 1, fc,. Note that Assumption 1 1 . 1 1 is verified whenever the kernels fi are bounded functions with 
support in a rectangle of the type B x ■ ■ ■ x B, fi{B) < oo. 



2 Discussion of the main results 
2.1 General bounds and mixed regimes 

Fix two integers d, m. Observe that, in the discussion to follow, one can take either d or m to be zero, 
and in this case every expression involving such an index is set equal to zero by convention. Our main 
results involve the following objects: 

- A vector = (Ai, A^) of strictly positive real numbers, as well as a random vector 

Xd = (X(i),...,x('*))^Po<i(Ai,...,A<i), 

that is, the elements of are independent and such that X^*^ has a Poisson distribution with 
parameter Ai, for every i = l,...,d. 

- A m X m covariance matrix C = {C{i, j) : i,j = 1, m}, and a vector N,„ — {N'^^\ iV'^'")) ~ 
^miO,C), that is, Nm is a m-dimensional centered Gaussian vector with covariance C. We will 
write H to indicate the {d + m)-dimensional random element 

i/ = (Xd,N„). (2.1) 

We shall also assume that X^ i N„i, where the symbol "X" indicates stochastic independence, 
and also that H JL rj, where rj is the underlying Poisson measure. 

- A vector — {F^^\ ....,F^'^^) of random variables with values in Z+ such that, for every i = 
1, d, i^W e domD and E(F,) = A,. 
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- A vector ~ (G^^-', G*^"'-') of centered elements of domZ). We use the notation 

F = (Fd,G™). (2.2) 
Note that, by definition, V is (T(?7)-nieasurable. 

Remark 2.1 Every asymptotic result stated in the present paper continues to hold if one allows the 
Poisson measure 77, as well as the underlying Borel measure space (Z, .2°, /i), to depend on the parameter 
n diverging to infinity. 

Our principal statement consists in an inequality allowing one to measure the distance between the laws 
of H and V . To do this, we shall need the following quantities, that are defined in terms of the Malliavin 
operators introduced above: 



ai(A<i,Fd) 

a2(F<i) 
a3(Fd) 



/3(Fd,G™) 
71 (C, G„i) 

72(G„i) 



I^i^fW (d^f(') - 1) D.L-^F^"^ fi(dz) 



i<i^j<d 



(df^'\-dl-^f^3^ 



fi{dz) 



+ E E 

d rn 

^^¥.l^\DL-^G^^^l\DF^'^ 
i=i ]=i 



D,F^^^D,F^''^D,L-^F^''> 



li{dz) 



= EiE 



c{j, k) - {dg^^\-dl-^g^^'^ 



/,(lh=1)^(ll-=-'='«l) 



li{dz). 



(2.3) 

(2.4) 
(2.5) 



(2.6) 
(2.7) 

(2.8) 



As we will illustrate in great detail below, the coefficients introduced in (I2.3p - (|2.8p should be interpreted 
as follows: (i) the sum ai(Ad, F,^) + a2(Fd) has the form X]i=i where each ai measures the distance 
between the laws of F^*) and (ii) a^{¥d) measures the independence between the elements of F^j, 

(iii) the sum 7i(C, G,„) +72(Gm) measures the distance between the laws of Gm and N,„, and (iv) 
/3(Fd, Gm) provides an estimate of how independent F^ and Gm are. Observe that A^ and C appear, 
respectively, only in ai and 71. Also, one should note the asymmetric roles played by Gm and F^ in 
the definition of /3(Fd,Gm). 
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Remark 2.2 A further connection between the quantity (j2.6p and the 'degree of independence' of 
and Grn can be obtained by combining the integration by parts formula of Lemma I A. 1 1 with the standard 
relation L = —SD, yielding that, for every j — 1, m and i = 1, d, 



E 



= E 



-dl-^g'^'\df^''> 



= Cov(G(j\i^W). 



A similar remark applies to the terms in a3{Fd). The fact that the dependence structure of the elements 
of the vector V can be assessed by means of a small number of parameters is a remarkable consequence of 
the use of the Stein and Chen-Stein methods, as well as of the integration by parts formulae of Malliavin 
calculus. In general, characterizing independence on the Poisson space is a very delicate (and mostly 
open) issue - see e.g. [^[501 [55] . 



We are now ready to state the main result of the paper, namely Theorem 12.11 The remarkable fact 
pointed out in its statement is that the above introduced coefficients can be linearly combined in order 
to measure the overall proximity of the laws of H and V. Observe that the estimate (|2.9I) involves an 
"adequate" distance di,{H,V) between the laws of the M''+™-valued random elements H and V. The 
exact definition of such a distance (which will be always a distance providing a stronger topology than 
the one of convergence in distribution on R'*+™) depends on the values of the integers d,m, as well as 
on the nature of the covariance matrix C, and will be formally provided in Section [3] (see, in particular. 

Definition 13.21 and Definition 13. 3p . For the rest of the paper, we will use the symbol " ''^ " to indicate 
convergence in distribution. 

Theorem 2.1 (Portmanteau inequality and mixed limits) Let the above assumptions and nota- 
tion prevail. 

1. For every d,m there exists an adequate distance d+(-, •), as well as a universal constant K (solely 
depending on \d and G ), such that 

d,{H, V)<K {ai(A<j, Yd) + a^i^d) + a^i^Pd) + /?(Fd, G„0 + 7i(C, G„) + 72(G„)} . (2.9) 

2. Assume Hn = (Fc(,n, Gm,n), n > 1, is a sequence of {d + m)- dimensional random vectors such 
that: (a) for every n, F^; „ is a vector of Jj^-valued elements of domD verifying 
Xi{n) := ¥,[Fn ^] — > Xi, (h) for every n, Gm,n — {Gn \ ...,G'"n^^) is a sequence of centered 

elements of domD verifying Cn{i,j) '■= ^[g[1^ Gn'^] — !• C{i,i) for i,i — l,...,m, and (c) as 
n — >■ 00, 

) + a2(Fd,„) + az{'Pd,n) + ^(Fd,„, Gm,„) + 7i(G„, G )^o, 

where Xd,n — (^i (n), Arf(n)), and C„ = {C„(i,j) : i,j = l,...,n}. Then, Hn V, where the 
convergence takes place in the sense of the distance di,{-,-). 

The proof of Theorem I2.1[ together with a detailed statement, is provided in Section [3?2] some direct 
applications of the mixed limit theorem appearing in Part 2 of its statement are described in Sections 
12.41 and |4l providing applications to random geometric graphs. Observe that the rest of our paper 
consists indeed in a series of applications of the estimate (|2.9p . obtained by properly selecting Ad, G, Fd 
and Grn- we will use this inequality to settle a number of open questions concerning probabilistic 
approximations on the Poisson space. The principal theoretical applications of Theorem 12.11 developed 
in the present work - namely to multidimensional Poisson approximations and asymptotic independence 
- are described in the next Sections |2. 2112. 31 
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Remark 2.3 Specializing (|2.9p to the case m = 1, d — 0, one obtains the main estimate in [37] . 
concerning normal approximations of Poisson functionals in dimension one. In the case m > 2, d = 0, 
(|2.9p coincides with the main inequality proved in |44| . where the authors studied multidimensional 
normal approximations on the Poisson space. Finally, the case d = 0, m = 1 corresponds to the one- 
dimensional Poisson approximation result proved in |35) . 

Remark 2.4 (About constants) Bv inspection of the forthcoming proof of Theorem [2.11 the constant 
K appearing in formula (|2.9p can be taken to be have the following structure: 

- If 771 = 1 and d > 1 (in this case, C is a strictly positive constant), 

1 - e-^- 1 - e-^^ 



^ 1 + 2V27r 
A = H h max 

C i=l,...,d 



A,,; 



A? 



where max0 by convention. 
If d> 1, 771 > 2, then 

K = 11 



max 
:=i,...,d 



1 - e-^' 1 



A. 



If d > 1 and 77i = 0, then 



K = 6 X ld>i + max 



1 - e 



-A; 



1 - e 



-A, 



A? 



(the case d = 1 follows from |35)). 



The values of the constants in the remaining cases (that is, when d is equal to zero) can be deduced 
form the main results of [371 HI] • 

We conclude this subsection with a refinement of Theorem 12. 11 -2 . providing useful sufficient conditions 
in order to have that the mixed term /3(Frf.„, Gm.n) converges to zero. 

Proposition 2.1 Assume H„ = Gm „), 77, > 1, is a sequence of {d + m)- dimensional random 



vectors such that Fd.n — {Fn^\ Fn""') is a vector of Ij^-valued elements of dom_D and G 

ice of centered elements of 
sufficient in order to have that lim /3(Fd_„, Gm,n) = 



{Gn \ ■■■,gI^'^) is a sequence of centered elements of domD. Then, the following two conditions are 



JziDM'^) l^{d. 



For every i = 1, d, the sequence n i— > 



There exists e > 1 such that, for every j = 1, 777, lim E 



is bounded; 



1+e 



li(dz) 



Proof. For every i,j, one can apply the Holder inequality to deduce that 



¥.{\DL-^G^^^\,\DF'^'^\) <E 



^{dz) 



X E 



D^L-'G^^ 



^{dz) 



(2.10) 



and use the fact that, since D^F^^'^ takes values in Z, then 



< 



for every e > 1. 
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2.2 Stable multidimensional Poisson approximations 



We will now discuss a class of multidimensional Poisson approximation results that are a direct conse- 
quence of Theorem 12. II Section r2.2.1l contains a general statement, whereas Section [2.2.21 will focus on 
sequences of vectors of perturbed multiple integrals. We will also establish several explicit connections 
with the multidimensional CLTs proved in |44j . 



2.2.1 General statements 

As indicated in the section title, with an additional small effort we will be able to establish limit theorems 
in the more general framework of stable convergence. The (classic) definition of stable convergence, in a 
form equivalent to the one originally given by Renyi in |52j (see also [2]), is provided below. 

Definition 2.1 (Stable convergence) Fix fc > 1. Let {X„} be a sequence of random variables with 
values in R''', all defined on the probability space {il, P) specified in Remark [1.2l Let X be a M'^-valued 
random variable defined on some extended probability space (r2',^',P'). We say that Xn converges 

st 

stably to X, written Xn — > X, if 



lim E 

n^oo 



E' 



(S) 



for every 7 £ R'^ and every bounded ^-measurable random variable Z. 



Remark 2.5 In this paper, we will be exclusively interested in stable convergence results where the 
limiting random variable X is independent of the cr-field This situation corresponds to the case where 
Z is defined on some auxihary probability space {A, £/, Q), and (fi', ^' , P') = {n x A,^ (g) £/,F (g) Q). 

Choosing Z = I'vci (S), one sees immediately that stable convergence implies convergence in distribution. 
For future reference, we now present a statement gathering together some useful results: in particular, 
it shows that stable convergence is an intermediate concept bridging convergence in distribution and 

convergence in probability. The reader is referred to [21' Chapter 4] for proofs and for an exhaustive 

p 

theoretical characterization of stable convergence. From now on, we will use the symbol — >■ to indicate 
convergence in probability with respect to P. 

Lemma 2.1 Let {Xn\ be a sequence of random variables with values in R'^. 

1. Xn ^ X if and only if {X^, Z) {X, Z), for every ^ -measurable random variable Z . 

st P 

2. If Xn — > X and X is ^ -measurable, then necessarily Xn — > X. 

3. If Xn X and {Yn} is another sequence of random elements, defined on (r2,^,P) and such that 
Yn ^ Y, then (X„,r„) 4 {X,Y). 

4- Xn ^ X if and only if (S) takes place for every Z belonging to a linear space % of bounded random 
variables such thatU^^'^^"^'^^ = L'^{VL,^,f). 

Remark 2.6 Properties such as Point 3 of Lemma 12 . 1 1 allow one to combine stably converging sequences 
with sequences converging in probability, and are one of the key tools in order to deduce limit theo- 
rems towards mixtures of probability distributions - e.g. mixtures of Gaussian random vectors. This 
last feature makes indeed stable convergence extremely useful for applications, for instance within the 
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framework of limit theorems for non-linear functionals of semimartingales, such as power variations, 
empirical covariances and other objects of statistical relevance. See the classic references 17[ and pTl 
Chapter 4], as well as the recent survey |46j. Outside a semimartingale framework, stable convergence 
on the Wiener space has been recently studied (among others) by Peccati and Tudor in [J^, Peccati 
and Taqqu [H], Nourdin and Nualart [3T] and Harnett and Nualart [18]. Some earlier general results 
about the stable convergence of non-linear functionals of random measures were obtained in [38, 39, 4(1,, 
by using a decoupling technique known as the 'principle of conditioning' - see \22 \ 163 ) . 



The next statement is a general stable multidimensional Poisson approximation result based on Theorem 
12.11 Recall that the total variation distance between the laws of two Z!j_-valued random elements A, B 
is given by 

dTv(A,B)= sup \F{Ae E) -F(B e E)\. (2.11) 

ECZ± 



A proof of Theorem 12.21 is detailed in Section 13.31 

Theorem 2.2 (Multidimensional stable Poisson approximations) Fix d > 1, let (Ai,...,Aci) S 



id let = (X(i), X((i)) ^ Po(Ai, Ad) be independent of rj. Let Fd.n = (Fn^\...,F, 



Ad) 



(■0 



n >1, be a sequence of -valued elements of domZ? such that E Fn 
^d,n = (Ai(n), ...,Xd{n)), n>l, and assume moreover that: 

) -|-a2(Fd,„) -l-a3(Fd,„) -> 0. 



A,(n) 



Xi . Write 



(2.12) 



Then, as oo, the law of Fd.n converges to the law ofK-d in the sense of the total variation distance, 
and relation \2. 9\) in the case m = provides an explicit estimate of the speed of convergence. If 
moreover. 



lim I 



D.Fi^^fiidz) 



= limE / D,Fi''>{D,Fi''> -1) fi{dz) 

J A 



lim E 

n—>oo 



[ D,L-^F,^:^fi{dz) 

J A 



= 0, 



(2.13) 



VI < « 7^ j < d, VA e jT^, lim E 



lim E 

n^oo 



fi{dz) 



(2.14) 



^{dz) — 0, 



then, Fd,n ^ 



Remark 2.7 1. Theorem 12.21 is the first multidimensional Poisson approximation result proved by 
means of Malliavin operators. In the case d = 1 (note that this implies — 0), the fact that 
condition (j2.12l) implies that (iry(Fi^„, Xi) ^ is a consequence of the main inequality proved 
in [35, . Applications of this one-dimensional result in random geometry appear in |35[ 159] . A new 
multidimensional Poisson approximation result in the context of random geometric graphs, based 
on the techniques developed in the present paper, appears in Theorem 12 .61 (c) . 
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2. A sufficient condition (that we will verify systematically in applications) in order to have that 
Oi2{Fd,n) + oi.3{Fd^n) ^ 0, is that the sequences 



(D,L~^F^^)\{dz) 



are bounded for every i and that, for every i ^ j, 



lim 



D,F^'\D,F^'^ ~l) fiidz)^ hniE 



H{dz) = 0. 



These conditions also imply that the middle term in (12.131) and the first term in (|2.14p are equal 
to zero. 

3. By a direct use of Point 4 of Lemma l2.1l (together with some adequate approximation argument), 
one can prove that another set of sufficient conditions in order to have stable convergence is that, 
for every A £ , every p > and every i = 1 , . . . , d. 



E 



fi(dz) 



0, 



where ...,Xp) = 1a{xi) ■ ■ ■ l^(a;p), and /q = 1 by convention. Albeit more easily stated 

than (|2.13p - (l2.14l) . these conditions are not simpler to verify in the applications developed in the 
present paper. 



2.2.2 The case of multiple integrals 

Now fix c? > 1. Our aim is to apply Theorem 12.21 in order to deduce a multidimensional Poisson 
approximation result for sequences of perturbed multiple integrals of the type: 

F,,„ = (fW,...,F('^)) = (xW +bW ),..., 4^^) +4'^) +/,,(/W)) , n > 1, (2.15) 

where : (i) each Fn"^ is a random variable with values in Z_|_, (ii) {xn : n > 1} is a sequence of positive 
real numbers, (iii) gi, . . . , (/d > 2 are integers independent of n, (iv) /^^ , . . . , /g^ indicate multiple Wiener- 
Ito integrals of respective orders qi,...,qd, with respect to the compensated measure fj, (v) for each 
1 < fc < d, fn''^ G Lg(/i*), and (vi) for each 1 < k < d, {B^^^ : n > 1} is a smooth vanishing perturbation, 
in the sense of the following definition. 

Definition 2.2 (Smooth vanishing perturbations) A sequence {Bn : n > 1} C L^(P) is called a 
smooth vanishing perturbation if i?„, L^^Bn G dom_D for every n > 1, and the following properties hold: 



lim E 


ll^-S„|||2(^) 


= hm E 


71— >00 


n— ^oo 


Um E 


ll^-S„||^4(^) 


= hm E 

n— J-oo 



\DL ^-Bn||i4(^) 



jim E[Bl] = (2.16) 

0, (2.17) 
0. (2.18) 

Note that, if (|2.17p - (|2.18l) are verified, an application of the Cauchy-Schwarz inequality yields that 



lim 



\DB„ 



= lim 



\DL-^Bn\\l: 



^ 
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Remark 2.8 Representing the Ornstein-Uhlenbeck semigroup as in [321 Lemma 6.8.1], one infers that 
the fohowing inequahties are always verified: 



> E 



\DL ^-B„|||2(^) 



> E 



\DL "^-Bn||L4(p) 



The following result is the announced multidimensional Poisson approximation result for perturbed 
multiple integrals. 

Theorem 2.3 (Poisson limit theorems on perturbed chaoses) Fix d > 1, \i, . . . , Xd > and let 

Xrf ~ Pod(Ai, . . . , Ad) be stochastically independent of rj. Define the sequence Fd^n, n > 1, according to 

Ii2.15\) . and assume that for each 1 < i < d, Xn^ — > Xi and E Iqiifn^)'^ — > Xi. Suppose also that: 

n—^QO L ^ J n— >-oo 

hm E[Fi^)i^,«] = lim {fi'\ fi'^) m,.^.) =0, I < i ^ j < d. (2.19) 

Assume moreover that the following Conditions 1- 3 hold: 

1. For every k — \, d, every r = 1, . . . ,qk and every / = 1, . . . , r A {qk — 1), one has that 

^ n—^oc 

2. For every k = I, d. the sequence n Wfn^'^WL^^fj^ik ) is bounded and, as n ^ oo, 



Z''k 



3. For every i ^ j such that qi = qj , 



lim ((/«)^(/(^■))^)^^^ =0. 



Then, 'Fd.n — ^ ~^d, o,rid the convergence of Frf_„ to takes place in the sense of the total variation 
distance. 



A proof of Theorem 12.31 is provided in Section 13.31 The following features of such a statement are 
noteworthy: 

- When specialized to the case d — 1, the assumptions of Theorem 12.31 coincide with those in [35l 
Theorem 4.1]. 

- In the case when qi ^ qj for every i ^ j, and apart from assumption (j2.19p . the statement of The- 
orem [53] does not involve any requirement on the joint distribution of the elements of the vectors 
Fd,n- This phenomenon mirrors some analogous findings concerning the normal approximation of 
vectors of multiple Wiener- Ito integrals on the Poisson space, as first proved in [il] . 

- In the case where qi — qj for i ^ j, Condition 3 in the statement follows automatically from 
(|2.19p . whenever fn^ and fn'' have the form of a multiple of an indicator function. 
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For the sake of completeness, in the next statement we present a slight refinement of the chaotic CLTs 
proved in (44j (the refinement resides in the stable convergence claim). Recall that the Wasserstein 
distance between the laws of two M™-valued random variables X, Y is given by 

dwiX,Y)= sup \E[g{X)]~E[giY)]\, (2.20) 
3eLip(i) 

where Lip(l) is the class of Lipschitz functions on R™ with Lipschitz constant < 1. 

Theorem 2.4 (Stable CLTs for multiple integrals) Fix m > 1, let N,„ = (A^^^), . . . , iV^")) - 
^(0,C), with 

C = {C{i,j) : i, j 1, . . . ,m} 
a m X m nonnegative definite matrix, and fix integers gi, . . . , > 1. For any n > I and i = 1, . . . , m, 
let gii^ e i^(/i'?*). Define the sequence Gm,n — {Gn \ ■ ■ ■ , Gn"''), 'n> I, as 

Gi:^^I,M^), n>l, z = l,...,m. 

Suppose that Assumvtion ll. 1\ is verified for every n, and also that 

hm E[GWgW] = X hm {gl:\ gl^^) = G(*,j), 1 < i,j < m. (2.21) 

Assume moreover that the following Conditions 1-2 hold for every k — 1, m: 

1. For every r — 1, . . . ,qk and every I — 1, ... ,r A (qk — 1), one has that 

2. Asn^ oo, /^,^ i^gi^^j d^l ^— ^ 0. 

st 

Then, 7i(G„,Gm,n) +72(Gm,n) 0, G,n^n N,„ and the convergence of G,n,n to N„i takes place in 
the sense of the Wasserstein distance. 

Remark 2.9 Apart from the covariance condition (|2.2ip . the assumptions appearing in the previous 
statement do not involve any requirement on the joint distribution of the components of the vector 
G„_m- Moreover, these assumptions are the same as those in [371 Theorem 5.1] (for the case m = 1) and 
[HI Theorem 5.8] (for the case m > 2). The somewhat remarkable (albeit easily checked) fact stated in 
Theorem 12.41 is that the same assumptions implying a CLT for multiple integrals systematically yield a 
stable convergence result. Note that this phenomenon represents the exact Poisson space counterpart of a 
finding by Peccati and Tudor [43] , concerning the stable convergence of vectors of multiple integrals with 
respect to a general Gaussian field. See [33j Chapter 6] for an exhaustive discussion of this phenomenon. 
CLTs on the Poisson space based on contraction operators have already been applied to a variety of 
frameworks - such as CLTs for linear and non- linear functional of Levy driven moving averages [ST] [40] , 
characterization of hazard rates in Bayesian survival models |12i 136] and limit theorems in stochastic 
geometry [13 dS]. 

We conclude this section by stating an application of Proposition 12. 1[ implying that vectors of (per- 
turbed) multiple integrals satisfying the assumptions of Theorem l2.3l and Theorem l2.4l are automatically 
independent in the limit. 

Proposition 2.2 Let the sequences {Fd^n ■ n > 1} and {G„i.n '■ n > 1}, respectively, satisfy the 
assumptions of Theorem \2.3\ and Theorem \2.4\ Then, l3{Fd.mGm^n) — >■ 0, as n ^ oo, and the two 
sequences are asymptotically independent. 

Several connected results involving [/-statistics are discussed in the next section. 
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2.3 Asymptotic independence of ^/-statistics 

We shall now apply the main findings of the paper in order to characterize the asymptotic independence 
of sequences of random variables having the form of U -statistics converging either to a Gaussian or a 
Poisson limit. Our basic message is that, under fairly general conditions, ?7-statistics verifying a CLT 
are necessarily asymptotically independent of any [/-statistic converging to Poisson. The criteria for 
Gaussian and Poisson convergence used below are taken from references [371 HH IS] and [55^ : to our 
knowledge, these references contain the most general conditions in order for a sequence of f/-statistics 
based on a Poisson measure to converge, respectively, to a Gaussian or a Poisson limit. 

By virtue of a de-poissonization argument borrowed from 116j . we will be able to deal both with pois- 
sonized and non-poissonized f/-statistics based on a i.i.d. sequence - see Proposition 12.31 The reader 
is referred to [26] for a survey of the classic theory of [/-statistics. See [3 19, 23, 60j, as well as the 
monograph [JS] and the references therein, for several examples of the use of [/-statistics in stochastic 
geometry. See [13l [27l [28l |35l [ST] [56l [57l |59] for new geometric applications based on Stein-Malliavin 
techniques. Albeit unified studies of Gaussian and Poisson limits for [/-statistics are available (see e.g. 
[23]), we could not find in the literature any systematic characterization of the asymptotic independence 
of [/-statistics in the spirit of the present section. 

Remark 2.10 Section contains another characterization of asymptotic independence of [/-statistics 
associated with random geometric graphs. Rather than using the general results discussed below, and 
due to the explicit nature of the kernels involved, we will establish such results by some direct analytical 
computations - allowing to obtain better rates of convergence, as well as results in higher dimensions. 

Since it is relevant for applications, we will explicitly work with a sequence of Poisson measures {ry„ : n > 
1}, each defined on the Borel space (Z, S") and controlled by a cr-finite measure /z„ ■possibly depending 
on n. Following [SU Section 3.1], we now introduce the concept of a U-statistic associated with the 
Poisson measure rjn. 

Definition 2.3 ([/-statistics) Fix k > 1. A random variable F is called a U-statistic of order k, based 
on the Poisson measure rjn with control /i„, if there exists a kernel h G L\{^^) such that 

F= '^W' (2.22) 

where the symbol 77^ ^ indicates the class of all /c-dimensional vectors x — (xi, . . . ,Xk) such that Xi £ rjn 
and Xi ^ Xj for every 1 < i ^ j < k. As formally explained in 51, Definition 3.1], the possibly infinite 
sum appearing in (I2.22|) must be regarded as the L^(P) limit of objects of the type J2xerj'' ^nAg 
q> I, where the sets Aq e Z'^ are such that n^lAg) < 00 and Aq Z'^ , a,s q 00. 

Example 2.1 (Poissonized [/-statistics) Assume {li : i > 1} is a sequence of i.i.d. random variables 
with values in Z and common non-atomic distribution p, and consider an independent Poisson random 
variable N{n) with parameter n > 1. Then, 77„(-) = X^t;^"^ '^^i (') * Poisson random measure with 
control fin = np- this framework, for every fc > 1 and any symmetric kernel h e Ll(pn) = Ll{(np)'^), 
the corresponding [/-statistic has the form 

F^ '^W= E h{Y,,,...,Y,J. (2.23) 

xgr;^ _^ l<il.,....ik<N{n);ii^ij 

The random variable obtained by replacing N(n) with the integer n in ()2.23|) is customarily called the 
de-poissonized version of F. 
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The following crucial fact is proved by Reitzner & Schulte in [511 Lemma 3.5 and Theorem 3.6]: 

Proposition 2.3 Consider a kernel h £ such that the corresponding U -statistic F in i2. 22]) is 

square-integrable. Then, h is necessarily square-integrahle, and F admits a chaotic decomposition of the 
form with 

/,(x,) /i,(x,) - j^^ ^ h{^,,^k-^)d^it\ X, e Z\ (2.24) 

for 1 < « < fc; and fi = for i > k. In particular, h = fk and the projection fi is in L].''^{ii\^) for each 
l<i<k. 

Remark 2.11 In [5T1 it is proved that the condition h E L^{iJ^) n L'^{iJ^-i) does not ensure, in general, 
that the associated f7-statistic F in (|2.22l) is a square-integrable random variable. 

The forthcoming Theorem 12.51 is the main result of the section. It is divided in three parts. Part 1 
collects some of the main results from |271 128j concerning the normal approximation of [/-statistics. 
Part 2 contains conditions for Poisson approximations of [/-statistics taken from [SSI Proposition 4.1]. 
Part 3 is new and states that, under the conditions appearing in the first two parts, any two [/-statistics 
converging, respectively, to a Gaussian and a Poisson limit are necessarily asymptotically independent. 

Remark 2.12 The bounds from [27l[28] stated below are easier to handle than the ones deduced in the 
seminal work [STj - albeit they are basically equivalent in several applications. The resulting conditions 
for asymptotic normality have been proved in |27| to be necessary and sufficient in many important 
instances. The conditions for Poisson approximations taken from [SO. should be compared with the 
classic findings of [23l |60] . 

Our framework is the following: 

- The sequence 

Gn ^ 5n(x), n>l, 

is composed of square-integrable [/-statistics of order k > 2 such that 5„ e L^{iJ^) n We 
write gi^m « = 1, k, for the ith kernel in the chaotic decomposition of (7„, as given in (|2.24p . We 
write (T^j = Var(G„) and write Gn = [G„ — E(G„)]/a-„. 

- For Gn as above, we set 

B{Gn;an) = ^ linax \\g,^n g^A + ll3^"l!L'l(p^) | > (2-25) 

where max ranges over all quadruples {i,j, r, I) such that l<l<r<i<j{i,j<k) and I ^ j 
(*) 

(in particular, quadruples such that I — r — i— j = ldo not appear in the argument of max). 

(*) 

- For an integer k' > 2, {An : n > 1} is a sequence of symmetric elements of 3f'' such that 

{An) < oo for every n. For every n, we define F„ to be the [/-statistic obtained from (|2.22p 
by taking /i(x) = hn{x) = fc'!^^l^^(x). To simplify the discussion, we may assume that each An 
is contained in a fc'-fold Cartesian product of the type Kn x • • • x Kn, with /i„(_ftr„) < oo, thus 
ensuring that each Fn is square-integrable. Accordingly, we denote by i = 1, k' , the ith 
kernel in the chaotic decomposition of Fn, and we also write A„ = k'\-^yJ^{An) = E[F„]. 
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- Define: 

p„ = sup^f,|(2/i, ...,yj) e : (yi, ...,j/j,ai, ...,ak'-j) G 
where the supremum runs over ah j = 1, fc' — 1 and all vectors (ai, ...,ak'~j) G Z'' ~^ .. 

Theorem 2.5 We denote by N and Xx, respectively, a ^{0, 1) and a Po(A) random variable, where 
A > 0. We assume that N ± Xx. 

1. There exists a constant Ck > 0, independent of n such that, 

dw{Gn,N) <CkB{Gn]<Jn). 

In particular, if B{Gn',crn) — > 0, then On converges in distribution to N, in the sense of the 
Wasserstein distance. 

2. There exists a constant Dk' > 0, independent of n, such that 

dTv{Fn,Xx) < |A„ - A| + Dk'^—^ (l + v/(A«+A2)(p„+p4) An. 

In particular, if An — t- 0, then Fn converges in distribution to Xx, in the sense of the total variation 
distance. 

3. Write Vn = {Fn,Gn), and H — {Xx,N). For an adequate distance, d^,, there exists a constant 
M = M{k; k'), independent of n such that 

d.{Vn, H)<Mx {a„ + B{Gn;cTn) + B(G„; 

In particular, if lim An = lim i?(G„;cr„) — 0, then Vn converges in distribution to H, and Fn 
and Gn are asymptotically independent. 

The next statement is the announced de-Poissonization result. 

Proposition 2.4 (De-poissonization) Let the notation of Theorem \2. 5\ prevail, and assume that, for 
every n, the Poisson measure rjn is defined as in Example \2.1\ Write F^^ and to indicate the de- 
poissonized versions of Fn and Gn- If lim An — lim _B(G„;it„) = 0, then G^) converges in 

n— foo n— >cxD 

distribution to H . 

2.4 Applications to random graphs 

We now demonstrate how to apply our main results to study multidimensional limit theorems for 
subgraph-counting statistics in the disk-graph model on R™. Our main contribution, stated in The- 
orem [32] below, is a new estimate providing both mixed limit theorems and multidimensional Poisson 
approximation results. The present section contains statements, examples and discussions; proofs are 
detailed in Section |4l Our notation has been chosen in order to loosely match the one adopted in [45l 
Chapter 3], as well as in [28l Section 3]. 

We fix m > 1, as well as a bounded and continuous probability density / on M™. We denote by 
Y — {Yi : I > 1} a sequence of R™-valued i.i.d. random variables, distributed according to the density 
/. For every n = 1, 2, we write N{n) to indicate a Poisson random variable with mean n, independent 
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of Y. It is a standard result that the random measure rjn — X^il;^"^ , where Sx indicates a Dirac mass 
at X, is a Poisson measure on M™ with control measure given by /i„((ix) = nf{x)dx (with dx indicating 
the Lebesgue measure on R™). We shall also write = — Mn, n > 1. Given positive sequences a„, 
we write ~ a„ to indicate that the ratio an/hn converges to 1, as n oo. 

Let {tn : n > 1} be a sequence of strictly decreasing positive numbers such that lim t„ = 0. For every 

n, the symbol G"(F;t„) indicates the undirected random disk graph obtained as follows: the vertices of 
G'{Y]tn) are given by the random set Vn = {Yi, 1^(71)} and two vertices Yi^Yj are connected by an 
edge if and only if \\Yi — yj||K'" G (0, t„). By convention, we set G'{Y\ i„) = whenever N{n) = 0. Now 
fix fc > 2, and let F be a connected graph of order k. For every n > 1, we shall denote by G'^iT) the 
number of induced subgraphs of G'{Y]tn) that arc isomorphic to F, that is: G'^{T) counts the number 
of subsets {ii, ik] C {1, N{n)} such that the restriction of G'{Y\ tn) to {Yi^ , Yi^} is isomorphic 
to F. Every graph F considered in the sequel is assumed to be feasible for every n: this means that 
the probability that the restriction of G'{Y;tn) to {Yi, Yfc} is isomorphic to F is strictly positive for 
every n. The study of the asymptotic behavior of the random variables G'^iT), as n goes to infinity, is 
one of the staples of the modern theory of geometric random graphs, and many results are known. The 
reader is referred to Penrose [45l Chapter 3] for a general discussion and for detailed proofs, and to [28l 
Section 3] and 51, Section 6] for some recent refinements. 

In what follows we will focus on the following setup: (i) ko,k are integers such that 2 < ko < k, (ii) 
the sequence {t„} introduced above is such that ~ n^"^^, (iii) Fq is a feasible connected graph of 
order ko, (iv) for some d > 1, (Fi,...,F(;) is a collection of non-isomorphic feasible connected graphs 
with order k. We also write 

G;(Fo)-E[GUFo)] 



Var(G;,(Fo))i/2 



The specificity of this framework is that, for such a sequence {tn}-, the random variables G^(Fo) and 
G'niTj) {j — verify, respectively, a CLT and a Poisson limit theorem. Our principal aim is 

to provide an exhaustive description of their joint asymptotic distribution. The following statement 
gathers together many results from the literature, mostly taken from |45[ Chapter 3] (for limit theorems, 
expectations and covariances) and [IQ (for the estimates on the Wasserstein distance). 

Proposition 2.5 Let the above notation and assumptions prevail. 

(a) There exist constants aq, &o > such that, as n ^ 00, 

E[G;(ro)] - ao«''"(C)''°"' ^ aon'^^-^"^'^''-^'^ ^ 00, 

anrf Var(G^(Fo)) bon'-''-'"'^/'-''-'^'^ 00. Moreover, the random variable G^(Fo) converges in 
distribution towards a o/K(0, 1) random variable, with an upper bound of order n^(''^''o)/'^{k-i) 
the Wasserstein distance. 

(b) There exist constants ai, > such that, for every j = 1, d 

E[GUF,)] ^ Var(G;(F,) ^ a,n'=(C)'-' ^ a,. 

Moreover, (G^ (Fi ),..., G'j(rrf)) converges in distribution to a d-dimensional vector (Xi, Xd) 
composed of independent random variables such that Xj has a Poisson distribution with parameter 
aj. 

(c) As n ^ 00, one has that, for every i,j ^ l,...,d, Cov(G;(Fo), G;(Fj)) ^ O („-('=-fco)/2(fc-i)) ^ 
and CoviG'niT,),G'niT,)) = O (n-VC^-D). 



18 



Remark 2.13 We could not find a proof of the multidimensional Poisson limit theorem stated at Point 
(b) of the previous statement. However, such a conclusion can be easily deduced e.g. from |45[ Theorem 
3.5], together with a standard poissonization argument. 

Plainly, Proposition 12.51 does not allow to deduce a characterization of the joint asymptotic distribution 
of the components of the vector 

Vn := (GUPi), GUr<i), GUPo)), n > 1. 

In particular, albeit Part (c) of such a statement implies that the random variables G'J^(ro) and G^(rj) 
are asymptotically uncorrelated for every j = nothing can be a priori inferred about their 

asymptotic independence. The following statement, which provides a highly non-trivial application of 
Theorem 12. 1[ yields an exhaustive characterization of the joint asymptotic behavior of the components 
of K. 

Theorem 2.6 (Mixed regimes in random graphs) For every n and every j = 1, ...,d, set Aj „ = 
E[G'„{Tj)]. Let iV ^(0, 1), Xd^„ = (Xi^„, ...,Xd,„) Pod(Ai,„, Ad,„), assume that N and Xd,n are 
stochastically independent, and write iJ„ = (X-d^m N). 

(a) There exist two constants A and B, independent ofn, such that, for some adequate distance di,, 

d*(F„, iJ„) < Ay/^ + Bn^^^ = O + n^ik^^ . (2.26) 

(b) Let Xrf ^ POc;(ai, ...,0^), where the constants aj have been defined in Part (h) of Proposition [^731 
be independent of N , and set H — (X.(i,N). Then, as n 00, Vn converges in distribution to H. 

(c) Write :— {Gj^(Ti), ...,G'„{Tci)), n>l. Then, there exists a constant C , independent ofn, such 
that 

dTvK,Xa,n) < Cy^^O [rr^^ . (2.27) 

Remark 2.14 (i) The estimates (|2.26|) - (j2.27|) and the content of Point (b) are new. We do not 
know of any other available technique allowing one to deduce the limit theorem at Point (b). 
Note that such a statement yields, in particular, the asymptotic independence of and G^(ro). 
The rate of convergence implied by formula (|2.27|) is probably suboptimal (the correct rate should 
be of the order of nt™ - compare with the statement of 051 Theorem 3.5] in the case of non- 
Poissonized graph). It is plausible that one could obtain a better rate by avoiding the use of 
the Cauchy-Schwarz inequality in the proof, and by estimating expectations by means of some 
generalized Palm-type computations (see e.g. [45l Section 1.7]). This approach requires several 
technical computations; to keep the length of the present paper within bounds, we plan to address 
this issue elsewhere. Previous classic references on geometric random graphs are [5| 123 1 [60 ] . 

(ii) A quick computation shows that if fc = fco -l- 1, the rate of convergence in (j2.26p is O <i('=-i) ^ 
and if fc > fcg + 2, the rate of convergence is O scs-i) ^ . 

Example 2.2 Let fco = 2, fc = 3, and consider the sequence of disk graphs with radius i„ such that 
t™ ~ n~'^/^. Define the following graphs: (i) Pq is the connected graph with two- vertices, (ii) Pi is the 
triangle and, (iii) T2 is the 3-path, that is, the connected graph with three vertices and two edges. Plainly, 
G^(ro) equals the number of edges in the disk graph, whereas G^(ri) and G^(r2) count, respectively. 
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the number of induced triangles and of induced 3-paths. Since Fi and r2 are non-isomorphic, Theorem 
12.61 can be apphed, and we deduce that G'j(ro), G'^{Ti) and G'„(r2) are asymptotically independent, 
and that they jointly converge towards a mixed Poisson/ Gaussian vector, with an upper bound on the 
speed of convergence of the order of n^^/^. 

We conclude this section by pointing out that an application of the de-poissonization Lemma 13.21 vields 
the following generalization of Theorem 12.61 The details of the proof are left to the reader. For every 
n, we denote by G{Y; tn) the de-poissonized random graph obtained as follows: the vertices of G{Y\ i„) 
are given by the random set Vn = {^i, Yn} and two vertices li, Yj are connected by an edge if and 

only if ||y,-r,|lR™ e (o,t„). 

Proposition 2.6 The conclusion of Theorem \2.6\ -fb) continues to hold whenever the disk graph G'{Y; tn) 
is replaced with the de-poissonized random graph G{Y;tn), and each counting statistic G'^iTi) is replaced 
by its de-poissonized counterpart. 



3 Proofs of the main theorems 
3.1 Preliminaries 

We will now introduce several classes of functions that will be used to define particular metrics used 
throughout the paper. We write g S C^(R™) if the function g : M™ — > M is bounded and admits 
continuous bounded partial derivatives up to the order k. Recall also the definition of the total variation 
distance drv given in p. lip . 



Definition 3.1 1. For every function g : 



R, let 

II II „ |gW ~9{y)\ 
Wghip ■■= sup— — , 

x^y \\x-y\\R^ 

where \\ ■ ||Rm is the usual Euclidian norm on R'". 
2. For a positive integer k and a function g G C^(M") , we set 



max sup 

l<ii<...<ifc<m x£R^ 



gk 



In particular, by specializing this definition to g^'^^ — g" and g^^^ — g'" , we obtain 

a2 



115 



max sup 



joo — max sup 

l<ii <i2 <i3 <m xeR"^ 



g3 



■9{x) 



dx^-^ dxi^ dxj 



-9{x) 



3. Lip(l) indicates the collection of all real-valued Lipschitz functions, from R to R. with Lipschitz 
constant less or equal to one. 

^. ^3 indicates the collection of all functions g G C^(R™) such that ||9||Lip 5: 1; lls'lloo 5: 1 O'f^d 

linioo < 1. 
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We now define the different metrics we will use. 



Definition 3.2 The metric between the laws of two Z'J^ x M- valued random vectors X and Y such 
that ¥,\\X\\2d E||y||^d < oo, written d^-^ (X, Y), is given by 

d^,{X,Y) = sup |E(/i(X)) -E(/i(y))|, 

where indicates the collection of all functions -0 : Z'J^ x R i— !■ M : (ji, . . . ,jd',x) i— )■ ■ ■ ■ ,jd',x) 

such that tp is bounded by 1 and, for all ji, . . . ,jd, the mapping x "0(^1, . . . ,jd',x) is in Lip(l). 

Definition 3.3 The metric between the laws of two Z'J_ x M'" - valued random vectors X and Y 
such that M\\X\\^d^^-jg^m, 'M.\\Y\\^d^y^m < oo, written d^^{X,Y), is given by 

d^,{X,Y) - sup \&{h{X))-nh{Y))\, 

where indicates the collection of all functions ip : Z^j. x R™ R : {ji, . . . ,jd;xi, . . . ,Xm) ^ 
■^Oii ■ • ■ : Jd; . . . , Xm) such that \ip\ is bounded by 1 and for all ji, . . . ,jd, the mapping {xi, . . . , Xm) 

tPUl, ■ ■ -JdlXi, . . .,Xrn) € %■ 

Remark 3.1 The indices 1 and 3 label the classes J^i and Jifs, respectively, according to the degree 
of smoothness of the corresponding test functions. The topology induced by any of the two distances 
, (i^3 is strictly stronger than the topology of convergence in distribution. 

We will sometimes need a useful multidimensional Taylor- type formula on Z"^. Given a function / on 
Z+, we write Af{k) = f{k + 1) — f{k), fc = 0, 1, and also A^/ = A(A/). More generally, given a 
function / : Z^ ^ R, for every i,j = 1, we write Aif{x^'^\ = , a;W + 1, - 

x^''')), and Af^ = Ai{Ajf). Of course, when rf = 1 one has that Ai = A and Af^ = A^. The 
proof of the forthcoming statement makes use of the following result, derived in ^35j Proof of Theorem 
3.1] (see also [4]). For every / : Z+ — > M, it holds that, for every k,a G Z+, 

f{k)-f{a)^Af{a){k~a)+R, (3.1) 

where i? is a residual quantity verifying 

|i?|<^^|(fc-a)(fc-a-l)|. (3.2) 

For the rest of the paper, we will use the following notation, which is meant to improve the readability 
of the proofs. If a; = [x^^\ . . . ,x^'^'') is an d-dimensional vector, for fc < p we will denote by x^'^''^^ the 
sub-vector composed of the fcth trough the pth component of x, i.e. x^*^'^^ — (x^'^^ . . . , x^^-*). Also, we 
set by convention x^^'^^^'^ = for every value of j . 

Lemma 3.1 Let / : Z^ ^ R. Then, for every x = (x^^' , a;^'')), a = {a^^\ ...,a('')) G U^, 

d 

fix) = f{a) + A./(a)(a;« - a«) + R, (3.3) 

i=l 
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where the residual quantity R verifies 



i<i^j<d 



Moreover, one has also the first order estimate: 

|/(x)-/(a)|< max ||A,/||ooX^| 

' ' 1=1 

Proof. Using p.ip . one has that 

d 

d 



where 



d 

\Ro\ < - . max ||A|/|U x ^ - a«||a:« - a« - 1|. 

Z 2— l,...,a ^ — ^ 



On the other hand, 

i—l i—1 i—1 

and formula (|3.3p is immediately obtained from the representation 

i=l 
z— 1 

as well as from the elementary inequality 

lAJia^^'^-'lx^^'"^) - AJ(a(i-^'),x(^+i-'^))| < \\Alf\\^ x \x^J^ - a«|. 
Formula p.4p follows from 

d d 



|/(2:)-/(a)| <^|/(a(i''-i),x('-'*))-/(a(i'^),x('+i''^))| < max ||AJ||o,^ 

— ^ ...a — ' 
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3.2 Complete statement and proof of the Portmanteau inequalities 

We provide below a precise statement of Theorem 12.11 including a discussion of the different cases in 
terms of dimensions and covariance matrices, each having its own associated metric. The technique of 
the proof is reminiscent of the computations contained in the classic paper :3j . For an explicit description 
of the constants K, Ki , see Remark 12.41 

Remark 3.2 We do not deal with the cases d = 1, m = and d = 0, m > 1 since they are already 
covered, respectively, by [35l Theorem 3.1], [37l Theorem 3.1] and [44l Theorem 3.3 and Theorem 4.2]. 
We could have dealt separately with the case where m > 2 and C > 0, by using a multidimensional 
version of Stein's method on the Poisson space, as done in [37l Section 3]: by doing so, we would have 
been able to consider test functions that are only twice differentiable, as well as bounding constants 
nicely depending on the operator norm of the matrices C and C^^. There is no additional difficulty in 
implementing this approach (albeit a considerable amount of additional notation should be introduced), 
and we have refrained to do so merely to keep the length of the paper within bounds. Finally, the results 
of [371 m] imply that, in the case d = 0, m > 1, one can drop the boundedness assumption for the test 
functions defining the distances djsri , 'i^s: well as the Lipschitz assumption for the functions compos- 
ing the class "^3. The forthcoming proof will reveal that these boundedness and Lipschitz properties are 
needed in order to deal with cross terms, that is, expectations involving both elements of F^ and Gm. 

Theorem 3.1 (Portmanteau inequalities: full statement) Let d,m be integers such that d\/m> 
1. Let H = (Xrf,N„i) and V = (Fd,Gm) be the (d + m) -dimensional random elements defined by 12.1]} 
and \2. 2\} respectively. Then, the following two statements hold: 

Case 1: d,m > 1. Consider the distance d^^ and d_^^, respectively, according as m = 1 or m > 2. For 
i — 1,3, there exists a universal positive constant Ki (having the form described in Remark \2.4\ l such 
that 

rf^. (V, H) < K, {ai{\d, Fd) + a2(Fd) + a^{¥d) + /3(Fd, G^) + 7i(C, G^) + 72(G^)} . (3.5) 



Case 2: d > 2, m — 0. In this case V = and H — X^, and one has that, for some universal 
positive constant K (having the form described in Remark \2.4\ l, 

drv (Fd, Xrf) < K {ai(Arf, F^) + a^iFd) + a^iF^)} . 

Proof. First of all, we observe that the conclusion of Case 2 follows from the computations leading to 
the proof of Case 1, by selecting a test function ip S .J^i uniquely depending on the first d variables. In 
what follows, K will denote a positive universal constant that may vary from line to line; by a careful 
bookkeeping of the forthcoming computations, one sees that such a constant K can be taken to have 
the form provided in Remark 12.41 

Now let ^p £ M'l. We want to deduce an upper bound for 

|E(^ (Frf,G„,))-E(V' (Xrf,N„0)|. 

We can assess such a quantity in the following way: 

|E(V'(Frf,G„))-E(7/.(Xd,N™))| (3.6) 
< |E (^ (Fd, G„)) - E (^ (Fd, N^)) I + |E (V- (F^, N^)) - E (^ (X.^, N„)) |. 
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The proof will consist of two main steps. In the first one, we will deal with E (tp (F^, Nm))— E (?/; (X^, N™)) 
and in the second one with E {ip (F^, G,„)) — E{ip (F^, N„i)). 

Step 1: Controlling the term E (-0 (F^, N,„)) — E ('0 (X^, N„i)). Such a term can be decomposed in 
the following way: 



E 



a 

(V (Fd, N„)) - E (V; (Xd, N„)) = J] E (x^i^^-D, Ff^''"^), N„) - V (x^'"), Ff^'+i^"^), N, 



fc=i 



We will now study separately each term appearing in the sum. In what follows, we write J^jj to indicate 
the probability measure given by the law of a given random element U\ integrals with respect to J^jj 
are implicitly realized over the set where U takes values. For any fixed 1 < fc < cZ, by exploiting 
independence, we have 



E 



E 



For a fixed {7S^'''-^\x^^+'^''^\y) G Z'f"^ x R™, we denote by 
the unique solution to the Chen-Stein equation 

verifying the boundary condition A^/(0) = 0, where :— (^z'''^'^"^\x^^\x'^^^^''^\y) . We recall 

(see e.g. [15]) that is given by /^(O) = /fc(l) - A/fc(2) and, for a; = 1, 2, 



fk{x) 



{x-l)\ 



E 



^(0H-E[V;(X«)] 



(3.7) 



(x-1)! 



Using the fact that 1-01 < 1 together with |T5j Theorem 2.3] and JTT] Theorem 1.3], we deduce that 
\fk\ < 3, |A/fe| < 2(1 - e-^'=)/Afe and |AV/c| < 4(1 - e-^'=)/A^0 Exploiting once again independence. 



^The upper bound on \fk\ can be reduced to 2 if one selects a solution of the Chen-Stein equation such that 
/(O) = 
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we can now write: 

e(v(x(1''=-i),f('=''^),N™) -e(v(x(1''=),f(W)^N„))) = (3. 



E [XuMk (x(i''^--i),F('=''^),N™j-(i?/fe(x(i^^-i),F(^-''^),N„,j 

Note that (since H is assumed to be independent of rf) in the previous expressions the Mahiavin operators 
act on random variables only through their dependence on the components of F^. We now need to 
explicitly calculate Dfk (X(1''=-i),F('='''),N„), and (by virtue of (HH)), one has 

D,/,(x(i^^-i),F('=^'^),N,„) 

- A- (x(i''=-i),F^'^),N„) - h (x(i-^-i),F('=^'^),N™) . (3.9) 

In order to deal with this quantity, one should first observe that, for every fc, the mapping /fc(-,N„i) : 
Z'' -> M given by 

a; H> /fc (x, N,„) 

takes values in [—3,3], and therefore \\Aifk{-,'^m)\\oo < 6 and || A| .j/fc(-, Nm)||oo < 12, for every i,j = 
1, d. One can now use Lemma l3.1l to deduce that 

DJ, (X(i^^-1),F(M),N„) = ^ A,A- (X(i''=-1),F('='^),N™) D^F'^^+Rf\ 

i—k 

where 

r . 

< 6 X <^ ^ \D,F^'^\\D,F^'^ - 1| + J2 \DzF^'^\\D,F'^^^\ 

I i—k ^^"^¥^3^^ 

Using the fact that (by definition) 

A,A (x(i^^-i),F('=^'^),N„) = AA (x(i^^-i),F('=''^),N™) , 

gathering the previous estimates together, and applying them to p.Sp finally gives: 

|E {4, (Fd, N,„)) - E (.A (Xd, N,„))| < {ai(Ad, F^) + a2(Fd) + a^{Ya)} . 

^; Controlling the term E (-0 (F^, G„j)) — E [ijj (F^, N„j)). This part is slightly more delicate, since 
one has to take into account the dependence between F^ and Gm- We have to consider two cases, 
namely m = 1 and m > 2. Note that, in the second case, it is not necessary to assume that the matrix 
C is positive definite. 
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(to = 1) In this case Gm and N™ are two real- valued random variables G E donil? and N ^ ^{0, 1). 
We will only consider the case C — 1, and one can recover the general statement by elementary consid- 
erations. For every x G 1J\_ and y G M, we write 

/^(a;,y) = e^'/2 T {^{x,a) - ^{i,{x,N)\}e-''^ da. 

J —oo 

It is well-known (see e.g. [331 Chapter 3]) that verifies the (parametrized) Stein equation 

dyU{x,y) - yUix,y) = ilj{x,y) - E[V'(x,iV)], yeR,xeZ'l, 

where we have used the symbol dy to indicate a partial derivative with respect to y. Moreover, thanks to 
the assumptions on ■0, one can prove that the following relations are in order for every x: \\f^{x, •)||oo 
V^, \\dyf^{x,-)\\ < 1, and \\dyyf^{x,-)\\oa < 2 (note that the partial derivatives dyyf^{x,-) are only 
defined up to a subset of M of measure 0). It follows that 

E iij (Fd, G)) - E (V- (Frf, N)) = E[dyU{Fd, G) - G/^(Fd, G)] 

= E[dyU{Fd, G)] - E[{-DL-^G, DU{Fd, G))l2(^^]. (3.10) 

Clearly, DJ^{Fd, G)^ A,^+B„ where 

Az '■— ftpiiFd)z,Gz) ~ /i>(F(i, Gz), Bz := ftjiiFd, Gz) — fip{Fd, G) 

Using a Taylor expansion as in ^STj Proof of Theorem 3.1], one sees that 

Bz=dyU{Fd,G)DzG + Rz, 

where |i?2| < {DzGY . Now observe that the mapping f^{-, Gz) : Z'' M : x i— )■ /^(x, Gz) is bounded by 
-\/27r, in such a way that || Ai/^(-, Gz)||oo < 2\/2tt, for every i — 1, ...,d. We can therefore use formula 
(|X4)) to infer that 

d 

\Az\<2V2^J2\^^^^'^\- 
Plugging these relations into p.lOp yields that 

|E {i, (Fd, G)) - E (7^ (Fd, TV)) I < if {/3(Fd, G) 71 (C, G) + 72(G)}. 



(to > 2) We use an interpolation technique analogous to the one appearing in [33] Proof of Theorem 
4.2]. For every t G [0, 1], we define 

$(0 E{V'(Fd, yi^G™ + ViN„,)}, 

in such a way that |E{?/;(Fd, G,„)} — E{-!/'(Fc;, N„i)}| < |<I>'(i)|rfi. Deriving with respect to t and then 
integrating by parts shows that 

^'[t)^At-Bt, 

where, with obvious notation. 
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and 



^ in 

Bt = -==^E[(-i^L-iG(j'),i^9j,^.V'(Fd,yr^G™ + V<N„))i2(^) 

.7 = 1 



2vr ^,=1 

where the random functions b^'^ and are given by 
and 

z ^ := 5j,^^(Fd, VT^t{Gm)z + VtN,„) - a;,^.7/;(F<j, VT^G™ + \/tN„). 
Reasoning exactly as in the proof of [44, Theorem 4.1], one proves that 



sup 

*6[0,1] 



At -Bp' <^{7i(^^,G„0+72(G„0}. 



1 



To conclude, we apply again Lemma l3.ll Start by observing that, since ^ 1 by assumption, one 

has that, for every i — 1, d, \\/S.idy.tp{-, VI — i(G„i)z + -\/tN„i)||oo < 2. We can now use (|3.4p to infer 
that 



i=l 



These estimates yield eventually that 



\Bf^\dt< 

10 V -L — I 

and the desired conclusion follows at once 



X /3(Fd,G„) =4/3(Fd,G™), 



3.3 Proof of Theorem [2721 



The first part of the statement is the same as Case 2 of Theorem 13. II In order to deduce the conclusion 
about stable convergence, one should fix an integer I > 1, as well as pairwise disjoint sets Ai, ...,Ai G i^, 
and then build an ancillary (d + Z)-dimensional vector 

Applying again Case 2 of Theorem 13. 1[ one proves immediately that conditions p.l2p - (|2.14l) imply that 
^'d+i n converges in distribution to (X^, 77(^1), f}(A;)). Since is independent of rj by definition, 
we deduce that, for every (71,..., 7^) e M'^, every collection Ai,...,Ai S ^ of disjoint sets and every 
random variable Z — ip{ri{Ai), ...,r/(A/)) with ip bounded, 



lim E 



= E[Z] X E 



An application of Point 4 of Lemma 12.11 yields the desired conclusion. 
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3.4 Proof of Theorem [2731 



Step 1: convergence in distribution. We start by proving that F^^n converges in distribution to X(j. 
Our plan is to apply Case 2 of Theorem 13.11 Exploiting the fact that each {Bn^} is a smooth van- 
ishing perturbation, and reasoning exactly as in the first part of the proof of [SS] Theorem 4.12], one 
sees that it is enough to prove that Conditions 1 and 2 imply that the five sums appearing in the 
definitions of ai (•), a2(-)j a3(') (see (|2.3p - (|2.5p ) all converge to zero, whenever one chooses the vector 
{iqiifn''), Iqdifn'^)) their argument. Again from the proof of [35l Theorem 4.12], we know that 
the assumptions in the statement imply that, for every i = I, d 



X.-q;'\\DI,^ifin'' 



lim <^ E 

Using the fact that the sequence 



{DJ,M^))\{dz) 



■E 



qfE 



{Dj,M^)r{Dj,M^)-ir^,{dz) 



^ 0. 



is bounded for every i, and by a standard application of the Cauchy-Schwarz inequality, we see that it 
is enough to prove that, for every i ^ j, 



lim <!E 

n—¥oo 



{DMfi^^)r{DJ,^{f^^^)rf^{dz) 



■E 



= 0. 



Using the computations contained in 37, p. 464], one sees that, for every i, 



2qi~2 



p=0 



(3.11) 



(3.12) 



where 



1 



fn \z, ■) *l f^'' (Z, ■){zi, Zp), 



9i — 1 r 

Gl^~^f^\z,-){zi,...,Zp) = ^ ^l{2g.-2-r-i=p}r! 

r=0 (=0 

(3.13) 

where the tilde indicates a symmetrization with respect to the variables represented by a dot (in such 
a way that the symmetrization does not involve the variable z), and the stochastic integrals are set 
equal to zero on the exceptional set composed of those z such that frP {z, •) ^J, (z, •) is not an element 
of ^^(/i^*'"^"'""') for some r,l. We can assume without loss of generality that qi < qj. Applying the 
isometric properties of multiple integrals using the Fubini theorem and integrating over Z, we see that 
the first summand in p. Ill) is a linear combination of objects of the type 

C„ = C„(/,r,s,t,p) := / {g\z,■)*ifr^:\z,■)Ji'\z,■)*tfi'\zr)}LH^.n^'(dz), 



where the indices verify the following constraints: (i) p = 0, 2^^ — 1, (ii) r — 0, gi — 1, (iii) I = 0, r, 
{iv)t = 0, ...,qj - 1, (v) = 0, ...,t, and (vi) 2q, - 2 ~ r - I ^ 2qj ~ 2 - t 

qi — qj^r — t, I — s — (and therefore p — qi — 1), one has that C„ — {[fn^)"^, {frPY^ 

In all other cases, one can prove that 



s — p. In the case where 
^ 0. 



\CJ < 



z)^0 
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(where the first inequality follows from the Cauchy-Schwarz inequality) by directly applying the com- 
putations contained in 37, p. 467], as well as the fact that (by assumption) sup ||/n*''||L'i(/j''i) < for 

n 

every i. We now focus on the second summand in (I3.1ip . We can use directly 011 Proposition 5.5] to 



deduce that, whenever qi — qj the quantity E {DIq.{fn^),DIq.(fn'))\2i^\ is equal to a finite linear 



combination of the squared inner product {fn \ fn'yj^2(^^qiy as well as of products of norms of the type 



qi—s{t,k) 



L2(^t + s(t.fc)) 



0')\2 



qj—s{t,k) 



J n 



(i)| 



(3.14) 



where s(t,k) = 2qi — k — t and the indices verify the constraints: k = l,...,2qi — 2, t = l,...,qi and 
1 < s(t, k) < t. On the other hand, when qi ^ q.j the same Proposition 5.5 in [44 implies that 



E 



DIq.{fn'),DIq.{fn ^))'j^2(^i^^) is & finite linear combination of products of norms of the type p.l4p . 
where s{t,k) — qi + qj — k ^ t and the indices verify the constraints: k — \qi — qj\, ...,qi + qj — 2, 
t = 1, qi A qj and 1 < s{t, k) < t. In both cases, the involved products of norms converge to zero 
whenever Condition 1 in the statement is verified, and Case 2 of Theorem l3 . 1 1 implies that Fd,n converges 
in distribution to in the sense of total variation. 



Step 2: stable convergence. We apply the second part of Theorem 12.21 In view of the previous compu- 
tations, and by reasoning again as at the beginning of the previous step, it is enough to show that, for 
every A € and every i — 1, d. 



This follows immediately from the relation 



J A 



where g{z) = 1a{z), as well as H/ri*' 5lli2(^,,-i) = (/« ' fn\g*l 5)l2(^2) (which follows from a 

Fubini argument). 



3.5 Proofs of Theorem 12.41 and Proposition 12.21 



Proof of Theorem \2.4\ In view of 44, Theorem 5.8] we only need to prove stable convergence. To do this 
we fix an integer d> 1, as well as disjoint sets Ai, ...,Aci E S'^. Using as before Point 4 of Lemma \TT\ 
the desired conclusion is achieved if we show that the (d + TO)-dimensional vectors (F^, G„j „), n > 1 
where = (77(^11), ri{Ad)), converge in distribution to (F^, N„i) (recall that N„i is independent of rj 
by definition). Define Ad = (^(^1), ^{Ad))- One has that ai{\d, F^) + a2(Fd) + a3(Fd) — 0, and also 
that, under the assumptions in the statement, 7i(C, Gm^n) + 72(G,„^„) (as a consequence of [44l 
Theorem 5.8]). To conclude, we have to show that (3{Fd,Gn.ni) — > 0. This follows immediately from 
Proposition 12.11 since the computations contained in ^37^ Proof of Theorem 5.1] imply that, under the 
assumptions in the statement. 



E 



^0, Vj-1, 



(3.15) 
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Proof of Proposition In view of Proposition 12.11 the conclusion is an immediate consequence of 
relation (|3TT5)) . ■ 



3.6 Proofs of Theorem 12.51 and Proposition 12.41 



Proof of Theorem \2.5\ For every n, let Xx„ be a one-dimcnsional Poisson random variable of parameter 
A„, and recall (see e.g. [1]) that dTv(^A,^A„) < |A — A„|. The distance in the statement can be 
chosen to be d,^^ (see Definition l3.2|) . An application of the triangular inequality and of the independence 
between Xx and yield that 

d4Vn,H)<dTv {Xx , Xx„ ) + d4Vn, (Xa„ , iV) ) . 
The conclusion follows from Theorem 13. 11 since one has that: 

- according to [SH Proof of Proposition 4.1], |A — A„| + q;i(A„, F„) + a2{Fn) < An] 

- according to (glj, 7i(l, G„) + 72(G'„) < CkB{Gn,(Jn); 

- by virtue of the Holder inequality, and of the fact that Fn takes values in , 



/?(F„,G„) <I 



\DA\%{dz) 



3 
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X E 






Hn{dz) 











<RX B{Gn;<7r 



Nl/2 



for some constant R independent of n, where we have used the fact that, since F„ and Gn 
both live in a finite sum of Wiener chaoses (see Proposition 12. 3p . then (a) the mapping n 

E \DzFnffinidz) is bounded, and (b) 



E 



DzL-^Gn 



Hn{dz) 



< 



DzGn 



^in{dz) 



< CkB{G„-(Tnf 



for every n. 



Proof of Proposition \2.4\ In view of the standard theory of Hoeffding decompositions (see e.g. [62]), for 
every n > 1, both F^ and (5° have the form of a [/-statistic of the type 



;=1 {ii,...,j,}C[n] 



(3.16) 



where [n] — {1, m is the order of the [/-statistic (so, m ~ k or m ~ k' , according as C/„ = G„ 

or Un — Fn), and every kernel [/„_/ is a symmetric function in I variables verifying the Hoeffding- type 
degeneracy condition: E[/7„^;(Yi, ...,Yi)\Yi, ...,Yi^i] — 0. The mean and variance of F„ and Gn are both 
converging, and this implies that, since the mapping 



71 = E[C/„]2 + ^ r )E[f/„,z(yi, Fz)"' 

1=1 ■ 



converges to a finite limit, then the sequences n (")E[[/„^;(Yi, Y;)^], I = l,...,m, are necessarily 
bounded. Now, it is easily seen that [/„ is the de-poissonized version of the poissonized [/-statistic [/^ 
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obtained by replacing [n] with [N{n)] = {1, iV(n)} in the second sum on the RHS of (I3.16p . The 
desired conclusion follows from the forthcoming Lemma |3.2[ whose proof uses computations from [161. 



Lemma 3.2 (De-poissonization Lemma) Let the above notation and assumptions prevail. Then, as 
n — )■ oo, 

Proof. Conditioning on N(n), and using standard results on the moments of Poisson random variables, 
yields (as n — ^ oo) 



/N{n) 

V I 



E[Un,i(Yi^ -.Yif] ^ c ~ lim E[U, 



1=1 

Conditioning again on N{n), we infer that 



n— >oo 



, -1 



where — X]^o ^ ''■^("^^)(") • conclude, it remains to apply the computations contained in 
[TBI p. 745], which imply that — )■ 1 for every I. ■ 

4 Random graphs: proof of Theorem 12.61 

The distance di, appearing in the statement is the distance d,^^ introduced in Definition 13.21 First of 
all we observe that, for every a ~ 0,1, ...,d, the random variable G^(ra) has the form of a LZ-statistic, 
that is: 

G'ni^a) = ^ ftr„,t„(a;i,...,xfcj, 

{xi,...,Xk^)eri'^'^^ 

where: (i) ka = k for a = I, d, (ii) the notation indicates that the sum runs over all ordered vectors 
{xi, ...,Xk^) such that each xi is in the support of ry„ and xi ^ Xf for I ^ I', and (iv) the quantity 
/ir„,t„ {xi, Xk^) is equal to l/Zca! if the restriction of G'(Y; i„) to {xi, Xk^} is isomorphic to Ta and 
equals otherwise. It is clear that, for every a, the mapping /ir„,t„ : (M™)'^'' R is symmetric and 
stationary, in the sense that it only depends on the norms \\xi — Xm\\R'", I 7^ m. We can now apply 
Proposition 12 . 31 to deduce that G[^{ra) admits the following chaotic decomposition 

G'^{^a)^E[G'^iTa)]+Y.I^{ha,^), (4.1) 
i=l 

where Ii indicates a multiple Wiener-Ito integral of order i with respect to the centered Poisson measure 
fin, E[G'„{Ta)] = /(R„)fc„ /Jr„,t„dMn% and, for i = l,...,ka, 

ha,n.i{xi, .■.,Xi) = i^'^] [ hY^.t„{xi,...Xi,yi,...,yka-i)l4t''^{dyi,...,dyka-i)- (4.2) 

Note that ha^n.ka — ^ra.t„- According to Theorem 12.11 our proof is concluded if we can show that 
the six quantities appearing in formulae (|2.3p - (|2.8p all converge to zero, as n — ^ oo, at a rate of 
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the order of O yn ^(fc-i) n nk-i) whenever one selects the following arguments: (1) = Vj[ = 

(G^ri), G'^{T,)), (2) G„ = Gi = G;(ro), (3) A, = A„,, = E[G;(r,)], I - 1, (4) C - 1, and (5) 
= fin- We already know from [SS] Sectfon 3] (see also Propositfon l2.5l -fa)) that the terms 7i(l, GJj(ro)) 
and 72(GJ^(ro)) both converge to zero at a rate r„ such that 



rn^O 



1 



Since n^°{f'^)^° ^ ^ n ^-i , this implies that we only have to focus on the remaining four terms. We 
start by analysing the term ai(A„, V^) and the first part of the term a^lV^). 



Select a,b — l,...,d. An application of the multiplication formula (|1.3p . together with the definition 
of the derivative operator and the representation (|4.ip . yields that 

{DGUra),~DL-'GUrb))m,^) 

= l{„=b}E[GUra)] 

k k lAj / . W { W ^ [ X\ 

Applying repeatedly the Cauchy-Schwarz inequality, one sees that, in order to prove that ai(A„,V^) 
and a-iiyl^) both converge to zero at the correct rate, it is sufhcient to show that, for every a, = 1, d 
and for every quadruple r, I) involved in the previous sum, 

\\K^n,^ A Vn,illL2(^j+J — = O [yf^^ = O (n^^T^) 

(the last equality is trivial). For any such {i,j,r,l) we define the function /i^ft'^'''' : (M™)" — >■ M, where 
a = a{i, j,r,l) — Ak^i— j — r + I, as follows: 

(^(r) ^(2) ^(3) ^(8)x, / (9) ^(6) (3) (S.). , . 

x/»r„,t„(Xj._.,x,_^,x^_j,X; )/ir,,t„(Xj._^.,x^-_^,x^_;,X; ), (4.4j 

where the bold letters represent multidimensional variables providing a lexicographic decomposition of 
(xi, Xq). For instance, one has that x^.'^^ = (xi, x^^^*^ = [xk-i+i, ■■■,Xk-r), and so on, in 
such a way that (x[,l^-, x|^'^, x|.5;, xj'^^ x[,^2^., xj*^'^, x[.'^-, x|*\ x[.^2^-) = (xi,...,Xq), and we set Xp"^ equal 
to the empty set whenever p — Q. Observe that each function h^^'j^'['^^ is bounded by l/k\^, and that the 
connectedness of the graphs ra,rb yields that the mapping (x2, ...,Xq) i— > /i^*'^'^''' (0, X2, ...,Xq), where 



stands for the origin, has compact support. Writing explicitly the squared contractions inside the 



integral, one sees that || /ia,n,i ^J. hb^njWi^i i+j-r-i-. is a multiple (with coefficient independent of n) of 



^ab't'^\^'i-' ■■■,Xa)f{xi) ■ ■ ■ f{Xa)dxi ■ ■ ■ dXa- 
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Applying the change of variables xi = x and Xi — tnUi + for i = 2, a, the above expression becomes 

a/,m\a-l / f/\ / ,(ij\r,/) 



n"(C)"' / / /ilTr(0,y2,...,ya)/(a; + t„y2)---/(x + t„y„)da;dy2---d2/a. 

JB-" J(R'")°-1 ' ' 

Since, by dominated convergence, the integral on the RHS in the previous equation converges to the 
constant 

f°'ix)dx / h^^f['^\{),y2,-.,ya)dy2---dya 
we deduce that \\ha.,n,i^i ^fc,nj ll^2(^.+i-'-') = O Since 

n (i„ ) =71 (i„ ) (ni„ ) 

and a — k > 1 for every possible choice of i,j,r,l, we immediately deduce the desired conclusion for 
cti{^n,Vn) ^^'^ the first part of 0:3(1^^). 

To deal with a2(Ki) and the second part of a3{V^), we apply the Cauchy-Schwarz inequality to write, 
for every a,b = 1, d, 



E / iD^G'jTa) {D^GU^a) - l)D,L-^G:,{n)\ fi^idz) < VA{a,n) x B{b,n), 

where A(a, n) = E i:>^G;(ra)2 (L>3G;(r<j) - 1)^ fi„{dz) and 

n) = E / [D,L-^GUrb)]''^ln{dz). 

One can easily verify that the sequence n 1— )■ B{b, n) is bounded (this is a consequence of the fact that 
G^Ff,) lives in a finite sum of Wiener chaoses). It follows that, in order to obtain the desired rate of 
convergence for this part, we just have to prove that, as n — ^ 00, A{a,n) = O {nt^). To do this, one 
applies again the multiplication formula (|1.3p (for every fixed z G R™) to deduce that, by virtue of the 
fact that hr^,t„ has the special form of an indicator multiplied by the factor kl^^, 

D,G'„ir,){D,G'^iTa) - 1) (4.5) 

z=l j=l r=0 ^ ^ ^ ^ (=0 ^ ^ 

^Ii+j-2-r-l{ha,n,i{z, •) *l ha,n,j{z, ■)) 
k-1 

-^t/t_i(/iQ,„,t(z,-)) := ^ C7(^;)- 

In the last equality, the set U represents the class of all indices (i, j, r, I) and t involved in the representa- 
tion of DzG'j^(ra){DzG'j^{Ta) — 1), whereas is the corresponding multiple integral process multiplied 
by the appropriate coefficient. To conclude, we apply the triangle inequality to deduce that 



E 

■yeu 



Sl{z)iiy,{dz) 



(4.7) 
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We will show how to deal with the quadruple (i, j, r, I) = {k,k,k — l,k — 1), which requires additional 
arguments than the others (which can be addressed in a straightforward way). In the particular case 
where (i, j, r, I) = {k,k,k — l,k — 1), we are looking at the term 



Thus, we have 



E 



^k,k,k-l±-l{z)lJ'n{dz) 



(z,-)] Unidz) 



kkl 



K,n,k{^^yi,-,yk-i)lJ'i ^{dyi,...,dyk-i)\ fin{dz). 



Using the fact that ha,n,k = ^r„,t„ along with the fact that hr^^t^ has the form of an indicator function 
and finally recalling the definition of /ia,n,i given by (j4.2l) . we can write 



E 



^k,k,k-l.k-li^)f^n{dz) 



(k - 1)! / hi„^,{z)^ln{dz) = {k- l)!|l/la.„,l|li^(p„) 



The analysis of the contraction carried out in the previous steps of the proof allow us to conclude that 
this quantity goes to zero at the correct rate when n goes to infinity (it corresponds to the (1,1,1,1)- 
contraction). Representing each remaining expectations in (14.71) as a contraction, and applying a change 
of variables analogous to the one described above gives the global and desired rate of convergence for 
^2(^1) ■well as for the second part of a3{V^). 

We now deal with the third and last part of Q!3(V^). Applying the same strategy, we can write, for 
every a,b^c — 1, ...,d with a ^ 



E 



\D,G'jra)D,GUn)D,L-'GUrc)\fin{dz) < ^Cia,b,n) x D{c,n), 



where C{a,b,n) = E J^^ D,G',,ir,fD,G'„{Tb)^finidz) and 



D{c, n) 



[D,L-^GUr,)]^^iJdz) 



Again, the sequence n 1— >■ D{c, n) is bounded and we can write, for a ^b, 



D,G'^{Ta)D,G'^{n) 

k k 



(4.8) 



1=1 



2—1 j—1 r—1 ^ ^ 



7e/ 



In the last equality, the set / represents the class of all indices {i,j,r,l) involved in the representation 
of DzG'^{Ta)DzG'^(rb), whereas C7 is the corresponding multiple integral process multiplied by the 
appropriate coefficient. This case is very similar to the previous one and the techniques used to prove 
that each of these expectation converge to zero as the correct rate are the same. However, there is 
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one additional term that was not present in the case of ck2(Ki)- This is the term corresponding to the 
quadruple r, I) — {k, k,k — 1, 0). We will detail this particular case. We have 



IE / Clk,k-i.oiz)^^n{dz) =E/ k^ll_^{ha,„,k{z,-)*l-ihb,n,k{z,-)) ^ln{dz) 

.JR" J JR"» 

= fc'*(fc-l)! / / hl^tJz,yi,...,yk^i)hl^t^^{z,yi,...,yk-i)fJ.^~^{dyi,...,dyk-^i)fin{dz). 

Jr™ J(R")'«-i 

Using the fact that /ir„,t„ and hr^^t„ have the form of indicator functions, we finally get 



E 



Clk,k^l.oiz)^in{dz) 



k^kl {h 



which is zero because Fq and Ff, are not isomorphic (^r„,t„ and hr^^t^ cannot be non-zero at the same 
time or Fq and Ff, would both be isomorphic to the same graph rendering them isomorphic to one 
another). This concludes the analysis of the term a3(V^'). 

It remains to deal with /3(V^, G'^(Fo)). Using relation (I2.10p with e = 3, we can write 



^(\DG'^{T,)l\DL-'G'^{T^)\) <E 



X E 



The term E 



is bounded and it remains to show that the term E 



DL-^G'„{To) 



goes 



to zero as n goes to infinity. For this, we will refer to |281 Section 3] where the rate of convergence of the 



term 72(G^(Fo)) is obtained by bounding it by a constant multiplied by yE DL ^G'^{To) ^n{dz). 

/ k-kp \ 

It is then showed that this last term goes to zero at a rate of O I n 2(^-1) j _ xhe difference here lies in 
the fact that the square root is replaced by a power j, yielding a rate of convergence of O I n ^'C'- 



When putting together all the rates of convergence for the different terms in the general bound, one sees 
that 



d^,{Vn, Hn) < A^nW^ + Bn ^'■^^-^) = O [n 2(^-1)+^ 
where A and B are positive constants that do not depend on n. This concludes the proof. 
Acknowledgments. We thank Christoph Thale for useful discussions. 
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A Appendix 



Throughout the Appendix, (Z, 3^) denotes a Borel space endowed with a non-atomic (T-finite Borel 
measure ^. We write r/ to indicate a Poisson measure on Z with control /i. As in the main text, 77 is 
assumed to be defined on some probabihty space (51, P) such that ^ is the P-completion of (j{r}). We 
also write ^^(p) = i2^fj,^,p). 

A.l Malliavin operators 

We now define some Malliavin-type operators associated with the Poisson measure 77. We follow the 
work by Nualart and Vives |33] . 

The derivative operator D. 

For every F E i^(P), the derivative of F, DF is defined as an element of ^^(P; L'^{fi)), that is, of the 
space of the jointly measurable random functions w : £7 x Z — > R such that E [J^ < 00. 

Definition A.l 1. The domain of the derivative operator D, written doniD, is the set of all random 
variables F e L'^{P) admitting a chaotic decomposition ([Ij such that 

fc>l 

2. For any F G domD, the random function z ^ D zF is defined by 

00 

D,F^Y.^h-i{fk(z,-))- 

k>l 

The divergence operator 5. 

Thanks to the chaotic representation property of 77, every random function u G L'^ {P , L'^ (fi)) admits a 
unique representation of the type 

00 

fc>0 

where the kernel fk is a function of fc + 1 variables, and fk{z, ■) is an element of L'l{^^). The divergence 
operator S(u) maps a random function u in its domain to an element of L^{P). 

Definition A. 2 1. The domain of the divergence operator, denoted by domS, is the collection of all 
u G L^{P,L^{ii)) having the above chaotic expansion (lA.ll) satisfied the condition: 

fe>0 

2. For u G dom(5, the random variable 6{u) is given by 

fe>0 

where fk is the canonical symmetrization of the fc + 1 variables function fk ■ 



39 



As made clear in the following statement, the operator S is indeed the adjoint operator of D. 
Lemma A.l (Integration by parts) For every G E dom_D and u £ domS, one has that 

E[GS{u)]^E[{DG,u)m^)]. 
The proof of Lemma [A. II is detailed e.g. in |34) . 

The Ornstein-Uhlenbeck generator L. 

Definition A. 3 1. The domain of the Ornstein-Uhlenbeck generator, denoted by domL, is the 
collection of all F G (P) whose chaotic representation verifies the condition: 

k>l 

2. The Ornstein-Uhlenbeck generator L acts on random variable F € domi as follows: 

LF--^fc/fe(/fc). 

k>l 

The pseudo-inverse of L. 

Definition A. 4 1. The domain of the pseudo- inverse of the Ornstein-Uhlenbeck generator, denoted 
by L~^, is the space Lq{V) of centered random variables in L^(P). 

2. ForF=Y. hifk) e Ll{V) , we set 

k>l 

L-'F=-Y,\lk{fk). 

k>l 

A. 2 Contractions 

Contraction operators play a crucial role in multiplication formulae and in the computation of expec- 
tations involving powers of functionals of the Poisson measure rj. In what follows, we shall define these 
operators and discuss some of their basic properties. The reader is referred e.g. to [42l Sections 6.2 and 
6.3] for further details. 

The kernel / ★J, g on ZP^'^^^^\ associated with functions / G LK^p) and g G where p,q > 1, 

r = 1, . . . ,p A q and Z = 1, . . . , r, is defined as follows: 

/ *r gill, Ir-h tl,,---, tp-r, Si,,..., Sg_r) (A.2) 

n'idzi, dzi)f{zi, , . . . , z;, 71, . . . , 7r-;, tl, , . . • , ip-r) 

^g{zi, , . . . , Zi,-fi, . . . Si, , . . . , Sq-r). 
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Roughly speaking, the star operator ' ' reduces the number of variables in the tensor product of / and 
g from p + q to p + q — r — I: this operation is realized by first identifying r variables in / and g, and 
then by integrating out I among them. To deal with the case I = for r = 0, . . . ,p A g, we set 

/ *r 5(T1i • ■ • J Trj ^Ij J • ■ • J tp—n Si, , • ■ • , Sq—r) 
= /(7l> ■ ■ ■ ,7r,^l, , ■ • ■,tp-r)9ill,- . . ,7r,Sl, , ■ ■ ■,Sg-r), 

and 

/ *0 9itl, ,---,tp,Sl,,...,Sq)=f(E) g{tl, ,...,tp,Sl,,...,Sq)= , . . . , tp)g{si, ,...,Sq). 

The kernel f -k^, g is called the contraction of index (r, /) between / and g. The above introduced 'star 
notation' is standard, and has been first used by Kabanov in (see also Surgailis ETj). Plainly, for 
some choice of /, g, r, I the contraction /^J, g may not be well-defined. The contractions of the following 
three types are well-defined (although possibly infinite) for every 1 < p < q and every pair of kernels 
g e i2(^P), / g L^^i): 

(a) f -kr 9(^1, --^Zp+q-r), where r = 0,....,p; 

(b) / 4 f{zi, Zq-i) = p{zi, Zq^i, for every / = 1, q; 

(c) for r = 0,....,p. 

In particular, a contraction of the type / -k^^ /, where I = 1, ...,q — 1 may equal -|-cxd at some point 
(zi, Zq-i). The following (elementary) statement ensures that any kernel of the type f^r 9 is square- 
integrable. 

Lemma A. 2 Let p,q > I, and let f G L'^ifJ-'^) and g € LKfi^). Fix r = 0,...,qAp. Then, f *r 9 ^ 

L-2^^^p+q~2ry 
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